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Abstract 

This paper studies rates of decay to equilibrium for the Becker-Doring equations with 
subcritical initial data. In particular, polynomial rates of decay are established when initial 
perturbations of equilibrium have polynomial moments. This is proved by using new dis¬ 
sipation estimates in polynomially weighted spaces, operator decomposition techniques 
from kinetic theory, and interpolation estimates from the study of travelling waves. 
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1 Introduction 


In this work we consider the Becker-Doring equations, namely the following (infinite) system of 
differential equations 




i = 2, 3,..., 

OO 

i=l 


( 1 . 1 ) 


where the Ji can be written as 


— UiCi(t)Cj(t) (t). 


( 1 . 2 ) 


and where (ad, [bi) are fixed, positive sequences, known as the coagulation and fragmentation 
coefficients respectively. 
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Polynomial decay to equilibrium for the Becker-Doring equations 


Becker-Ddring systems form a subclass of the more general coagulation-fragmentation equa¬ 
tions. In typical physical applications the Ci represent the discrete distribution function of par¬ 
ticles of size i, and the evolution given by (HU represents the mean field approximation of the 
evolution of the distribution function q. In particular, Ji{t) represents the net rate that parti¬ 
cles of size i and size I either join to form particles of size i -I- 1, or conversely are emitted by 
spontaneous breakup. Thus we are primarily interested in positive solutions, whose first moment 
is preserved in time, meaning that 

CO 

Ci > 0, '^ici{t) = g{t) = g, Vt > 0. (1.3) 

i=l 

The Becker-Doring equations are used to model reactions in various physical settings, such as 
vapor condensation, phase separation in alloys and crystallization. This model was first proposed 
in [3] , and was modified to the form we are considering in [5],[in]. A good mathematically-oriented 
review can be found in [18) . 

The well-posedness and convergence properties of the Becker-Doring equations have been 
well-studied. In particular. Ball, Carr and Penrose [3] demonstrated the existence of “mass”- 
preserving, non-negative solutions to this system, namely solutions of (HU satisfying (HU- A 
later work |13j established well-posedness (including uniqueness) for any initial data with finite 
first moment, namely the space where the “mass” is well-defined. Ball et al. [2] also demon¬ 
strated that as t —?> oo solutions must converge to some equilibrium (Qi), where (Qi) is uniquely 
determined by g. Furthermore, they prove the existence of a value Ps such that if p < Ps then 
the convergence to (Qi) is strong. On the other hand, if p > Ps then there is a loss of mass to oo, 
and the convergence is only weak. Any initial data satisfying p < ps is called subcritical, while 
data satisfying p > ps is supercritical. 

The goal of this paper is to quantify the trend to equilibrium in the subcritical case (p < Ps). 
Specifically, our goal is to establish uniform, local rates of convergence to equilibrium in spaces 
with polynomial moments. 

We define the detailed balance coefficients, a sequence {Qi), by the equations 

Ql 1 ; Qiei'i , 'I (^■^) 

We note that the equilibrium solution Qi can be written as 

Q^ = Q,z\ (1.5) 

where the parameter z is related to the mass p in the subcritical regime through the equation 

OO 

^^iQi = p. 

We note that Qs is linked to the radius of convergence Zs of the power series with coefficients Qi. 
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Part of our interest in studying these equations is precisely that we believe that the Becker- 
Doring equations are a suitable prototype of more general coagulation-fragmentation equations 
with detailed balance. Indeed, we suspect that many of the interesting phenomenon that occur 
for the Becker-Doring equations may be typical of other systems with detailed balance. 

Convergence to equilibrium was proven by Ball, Carr and Penrose [5] using an entropy func¬ 
tional. Specifically, they prove that the quantity 

OO / \ 

V'(c):=gc.(l„g|--l) 

is weak-* continuous and that V{c{t)) is strictly decreasing. 

Later, Jabin and Niethammer m proved an entropy dissipation inequality which gives a 
uniform dissipation rate for regular data. In particular, they proved that if the initial data 
decays exponentially fast, then the solution converges to equilibrium with a rate bounded by 
in the mass-weighted space. 

In a recent work, Cahizo and Lods [6] improved this bound to . They do so by observing 
that the Becker-Doring equations da have a type of symmetric structure, which we make use 
of below. In particular, if we write the Becker-Doring equations in terms of a perturbation of the 
equilibrium solution 

Ci=Qi{l + hi), (1.6) 

then we may express the original equation m in the form 

jh = F{h^{t))h, 


and the mass constraint (11.31) as 

OO 

0 — ^ ^ Qi'i'hi. 

2=1 

(1.7) 

We note that 


F{g) = L + gP, 

(1.8) 


where L and P are both linear operators. Cafiizo and Lods rewrote the operator F{g) in weak 
form, satisfying 

GO OO 

2^1 2^1 

for all in a suitable space of test sequences. They then linearized the equation and consider 
the operator L, which is given in weak form by 

OO OO 

^ ^ Qi (T/i)202 ^ ^ CliQiQli^hl T hi ^2-1-1 ) (02+1 02 01 ) ■ 

2=1 2=1 


(1.9) 
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If we consider an P' space weighted by Qi then this form is clearly symmetric. Additionally, if 
(ci) is a solution of (HU and (hi) is determined by (11.61) then we have that hi € [—l,oo) and 
that Qiihi = 0. We thus define the Hilbert space H by 

H ■.= \{hi) :\\h\\pi^Q.) ■= < oo, Qjihi = 0 


with the natural induced norm || • \\h = || • \\p{Qi) and inner product (•,•)//■ Cahizo and Lods 
demonstrated that the linear part (L) of the Becker-Doring equations has a good spectral gap in 
H, or precisely that for some constant Ac > 0 the following holds, independent of h: 


Lh'^ ^ ^ OiiQiQ\{h\ hi ^ Ac (h, h^j^ . 


( 1 . 10 ) 


i=l 


A key point is that the mass constraint (HU precludes the null vector hi = i. Detailed quanti¬ 
tative estimates of Ac can then be obtained using Hardy’s inequality—see [6] for details. 

Cahizo and Lods then utilized a priori bounds from HD to control the non-linear term and 
establish a rate of convergence to equilibrium. More precisely, defining the Banach space 

:= |(hi) : ||h||^i(Q.er,i) := < cx), = 0 < 77 < 1, 

with the induced norm || • ||x, = || • they prove that for 0 < 77 < 77 , given initial data in 

Xfj then the solution must converge at a uniform exponential rate in . A key technical aspect 
of their proof was an operator decomposition technique from [9], which permits an extension 
of the spectral gap of L from H to A^. We recall (see [6]) that the space H is continuously 
embedded in A^ for 77 > 0 sufficiently small, precisely because the Qi are exponentially decaying. 

Our aim in the present paper is to analyze the trend to equilibrium for a wider class of initial 
data, for which the a priori bounds from m are not available. We dehne the Banach spaces 

A/i, := ^ (hj) : := ^ ) Qji \hi\ K. 00 , ^ ) Qjihj = 0^ , ^ ^ 1, 

with norm || • = || • ||fi(Q.i»=). The main result of our paper is as follows. 

Theorem 1.1. Let {hi(t)) defined by (11.61) represent the deviation from equilibrium of a solution 
{ci(t)) to the Becker-Doring equations (see Definition \l.‘A) . Assume that the model coefficients in 
(HU satisfy I|1.11I) - (I1.14I) below. Let m and k be real numbers satisfying m > 0 and k > m -\- 2. 
Then there exists positive constants Sk,m,Ck,m so that if ||h(0)||js:j^^ < 6k,m then we have that 

\\h{t)\\x,+^<Ck,m{l + t)-^'‘-"^-^'>\\hm\x,+, forallt>0. 


In order to obtain this result, we establish detailed estimates on the semigroup generated by 
L in the spaces Xk by using new dissipation estimates, together with the spectral gap estimate 
(11.101) . the operator decomposition result from [9] and interpolation techniques from Engler’s 
work on travelling wave stability [ 8 ]. We then prove a local stability result in Xk (see Theorem 
EU), which along with Duhamel’s formula proves the desired result. 
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1.1 Assumptions and Preliminaries 

We impose the following assumptions on our model coefficients: 


Qi 

> Cl > 0 Vi > 1, 

(1.11) 

y Oi + l 

lim - 

= 1, 

(1.12) 

2—)-oo 

=: — G (0,oo) 

T ai 

lim — 

(1.13) 

2—^00 bi 


Qi , bi 

< (72* V* > 1, 

(1.14) 


with Qi , bi as in (HH) and and where Ci, C 2 are fixed constants, independent of i. 

Following [2], we define a solution to the Becker-Ddring equations in the following way 

Definition 1.2. A function {ci{t)) is a solution to the Becker-Ddring equations on [0,r) if 

1. *|ci| < 00 for all t G [0,T). 

2. For all i we have that Ci{t) is continuous in time, and non-negative. 

3. The following (well-defined) equations are satisfied 

Ci(i) = Cj(0) + / {Ji-i{s) - Ji{s)) ds, i>2, 

Jo 

Ji(s)^ ds. 

Throughout the paper we will be considering solutions {ci(t)) of the Becker-Doring equations 
(ED with some fixed, subcritical mass, meaning that for some 2 : < Zs, we have that the Qi 
defined by (11.51) will satisfy 

00 00 

'^QA = g = 

i=l i=l 

Using p.4l) . (11.51) and (11.131) . it is immediate that 

lim = _L < 1 , (1-15) 

i-^oo Qi Zs 

This naturally implies that the Qi are exponentially decaying. 

Also, by combining (I1.12|) and (jl.l3|) . we observe that 

ai{z-\-6) = ai{Qi-\-6) < bi, yi>Nz, (1-16) 

for some (5 > 0 and that are fixed and independent of i, but possibly dependent on 2 . All 
of these assumptions are fairly standard, and versions of them can be found in [TJ [SJ 111) . In 


ci(t) = ci(0) - / (ji(s) + ^ 
•^0 V i=i 
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particular we note the similarity of (11.161) with the assumption given in [T]. In that work Ball 
and Carr make the assumption that 

a-iZ < bi 

for i > N, and for all z < Zg. In that work, this assumption was made in order to guarantee 
that V{c(tn)) converges to the minimum value of V, where I^ is a suitable entropy functional. 
In their work, coefficients were required to be 0{i/ log(i)), but this was subsequently relaxed in 
m- These assumptions were also utilized in and m 

One of the primary advantages to our method is that it lays bare the mechanism causing con¬ 
vergence to equilibrium. Inequality (jl.l6p arises naturally in attempting to establish dissipation 
estimates, thus motivating the analytical need for such assumptions. More importantly, (11.161) is 
satisfied by many of the relevant physical models. For example, one physically-motivated form 
of the model coefficients is (see [16]) 

= bi = ai (^Zg +, ae(0,l], /xG[0,1], q > 0. 

For this model we have 

bi QlCii ^ (^s 

which naturally implies that assumption (11.161) is only satisfied in the subcritical setting. 

2 Linearized stability estimates In Xi 

In this section we establish stability estimates for the semigroup generated by the operator L, 
in the space Xi. Following , our goal is to use an operator decomposition technique to derive 
uniform bounds on in Xi. This technique was first developed by Gualdani, Mischler and 
Mouhot |9] to study the Boltzmann equation, and was previously applied to the Becker-Doring 
equations by Canizo and Lods |5]. Here we generalize the technique to the case of evolution 
families. 

We remark that the symbols M and C, with various subscripts, will represent generic con¬ 
stants whose value may sometimes vary line by line. The notation C{Y^ Z) denotes the space of 
bounded linear operators from Y to Z, and C{Y) = C{Y, Y). In this paper, the term “semigroup” 
always refers to a strongly continuous semigroup of linear operators. 

Proposition 2.1 (Extension Principle). Let Z G Y be Banach spaces, with Z continuously 
embedded into Y. Let I = [0,r) with T = oo permitted, and let {H(t)}tg/ and {B{t)}t^i be 
families of linear operators on Y. Suppose that 

1. {A{t) + B{t)}tei generates an evolution family on Z, satisfying 

\\U^{t,s)\\ciz) < forO<s<t<T, 


a for some G K. 
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2. B{t) is “regularizing,” meaning that B{-) € C{I]C{Y,Z)), and that ||S(t)||£(y,z) < Mb, 
uniformly for t G I. 

3. generates an evolution family V on Y, satisfying 

ll^(i,s)ll£(y) < for 0 < s <t <T, 


with Xz > Xy- 

Then {A(<) + B{t)}t^j generates an evolution family on Y with bound 


\\U^it,s)\\ciY) < forO<s<t<T. 


( 2 . 1 ) 


Proof This proof is almost identical to that of Theorem 3.1 in [B], with the necessary changes 
to the setting of evolution families. We give the proof for clarity. 

First, it is evident that {A{t) + B{t)}t^i generates an evolution family since B{t) is bounded 
and continuous in t (see, for example. Theorem 5.2.3 in [14)1. Thus our goal is to prove (EH). 

Using Duhamel’s formula we can write the evolution family generated by A{t) + B{t) as 
follows: ^ 

{t,s)h{s) = V{t,s)h{s) + f {t,r){B{r)V{r,s)h{s)) dr 

J S 

This formula can be rigorously justified in the current setting by applying Lemma 5.4.5 in 
M- We then estimate 


||C/^(t, s)/i(s)||y < Mye 


— \Y(t — s) I 


ll^(s)lk+ / \\U^ {t,r)B{r)V{r,s)h{s)\\Y dr. 


As B maps from Y to Z we can replace with inside the integral, and then estimate 
using the decay estimate in Z to infer 


\\Uit,syh{s)\\Y <Mve 


-Ar (t-s) 


||/i(s)||y+ / Mze 


-\z{t-r) 


\\B(r)V{r,s)h{s)\\z dr. 


Using our bounds on B and V we obtain 


\\U^{t,s)h{s)\\Y < Mve-^^^^-^^\\h{s)\\Y + ms)\\YMvMzMBe-^^^^-^^ 
<MYe-^-(‘-^)||h(s)||Y, 


g-(Az-Av^)(t-r) 


which is the desired result. 


□ 


Remark 2.2. When A and B are constant in time this reduces to a statement about semigroups, 
and indeed in that case the statement and proof are found in [Q . In this section we only use the 
proposition to prove bounds on the semigroup but in Section |4] we will use it in the case of 
evolution families. 
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We emphasize that the previous result is valid when Ay = 0, meaning that the result applies 
to semigroups that are only stable. 

Next, recall that the operator L is determined by the weak form (11.91) . We write 

L = A + B, 

with the operator A determined via the weak form 

OO CX) 

^ ^ ^ ^ QiQlOji{h^ ) (02y 1 (pi Q N —iQ l^N ^l); 

i=l i=N 

( 2 . 2 ) 

where we fix some fV > fVz + 1, with given in (11.161) . We take the domain of definition for 
both A and L initially to be the set of sequences with finite support that satisfy (IE3, namely 
having zero “mass”. We note that if we set (j>i = i we get zero, implying that A and B both map 
into the space of sequences with zero mass. 

We first give an elementary bound on L and P, which indicates a minimal size for the domain 
of the closure of these operators. We will subsequently show that B is bounded, which in turn 
means that this also gives information about the domain of the closure of A. 

Lemma 2.3. For any m > 0, and for some constant Cm the following bound holds 
||r/l||xi+„ < Cm\\h\\x2+m \\Lh\\xi+^ < Cm\\h\\x2+r^- 

Proof. We only show the estimate for L, as the estimate for P is essentially identical. We simply 
estimate 

C30 

\\Lh\\x,^^=Y.Q^{Lh)A^+'^^sgn{Lh), 

i^Q 

OO OO 

< ^ ^ Qi{ciiQi bi)\hi\3{i 1 )^+^ _|_ 3 |/^^| 'y ^ QiQiai{i 

i=l i^l 

OO 

i=l 

where we have used (|1.14l) . This proves the lemma. □ 

In order to use the extension principle, Pror)osition l2.11 we first prove that B is “regularizing.” 
(Recall H cXi.) 

Lemma 2.4. The operator B is a bounded operator from Xi to H. 

Proof. We compute in weak form: 

OO N — 2 OO 

2=1 2=1 2=1 

+ QN-iQiax-ihx-iif'N — (fx-i — 4>i) 

=: Bi{h,4>) + B2{h,4>) + B3{h,cj)). 
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By Cauchy-Schwarz, the fact that 0 < c < QijQi+i < C < oo by (I1.15L and the equivalence of 
finite dimensional norms, 


\B,{h,4’)\<C 


/N-1 


1/2 




^Af-l 


- 1/2 

Q^hA <c\mH\\h\\x,. 


Furthermore, 


1/2 


1/2 


\B2{h,(^)\<C\h,\[Y,Q^aU 






Similarly, |_B 3 (li, (^)| < CH/iUxi ||<^||_f/- Setting 0 = Bh then establishes the desired result. 


□ 


Next we need to show that A, or more precisely its closure, generates a contraction semigroup 
on Xi. We recall the following definition from Pazy [14]. 

Definition 2.5. Let x G X, with X a Banach space. Define 

J{x) := {a:* e X* : {x*,x)^, ,, = ||x||^ = ||x*|||.} . (2.3) 

A linear operator A with domain of definition dom(A) G X is called dissipative if for every 
X G dom(A) there exists an x* G J(x) such that 

{x , X* ,x — ® 

By way of notation, when X = ii{QiWi) and ||fi||x = X^fci Qi'>J^i\hi\ we will write 

OO 

{sgn{h),(l))x,^x ■= X! QiWi(l>iSgn{h^). 

i=l 

By the definition of J{x), namely (12.31) . it is clear that if {sgnih), Ah) x — ^ ^ 

domain of definition of A then A is dissipative. 

Proposition 2.6. The operator A given by (12.21) is dissipative on Xi. 

Proof. Rearranging our sum and using dn to say QiQiOi = Qi+ibt+i, we find that 


(sgn(/i),A/i)x*,Xi 

OO 

= ^ QiQiaihi{{i -\- l)sgn{hi^i) -isgn{hi) -sgn(/ii))) 

i^N 

oo 

- X! Qihhi{isgn{hi) - {i- l)sgii(/ii_i) - sgn(/ii)) 

i^N 

oo 

= ^ Qihi(^Qiai{i + l)(sgn(/zi+i) - sgn{hi)) + bi{i - l)(sgn(/zi_i) - sgn(/^i))^ 

i^N 

oo oo 

“h ^ ^ Qi\hi\{(^iQl H“ Sgll(/ii) ^ ^ Qihiipi QlO/i) 
i^N i=N 


—: El H- E2 + -E’s, 
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Because hi{sgn{hi±i) — sgn(hi)) < 0, we see Ei < 0. Furthermore, we have that 

OO 

E 2 + E^ = 2 'y ^ Qi\hi\{aiQi — bi) < 0 . 
i=N 

sgn(hi)^sgn(hi) 

This readily implies that A is dissipative (see Definition 12.5|) . □ 

Remark 2.7. In the case that ~ i it is probably possible to prove that L has a spectral gap 
in Xi. We do not pursue this line of analysis, because in most of the physical cases ^ i“, with 
a G (0,1). 

Next we recall two results from [5] (Corollary 2.11 and Theorem 3.5), that concern the closure 
of L (which we also denote below by L). 

Proposition 2.8. For some > 0, the operator L generates a contraction semigroup on El 
satisfying 

\\e^*\\c{H) < e for all t >0. 

Furthermore, for 77 > 0 sufficiently small there exists constants M and Xrj > 0 so that the operator 
L generates a semigroup on satisfying 

\\e^*\\cix„) < Me~^A forallt>0. 

Next we prove that the closure of A indeed generates a semigroup. 

Lemma 2.9. The closure of A (which we also denote by A), generates a contraction semigroup 
on Xi. 

Proof. The Lumer-Phillips theorem (see e.g. [7] Theorem II.3.15) states that the following are 
equivalent for a densely-defined, dissipative operator A: 

1. The range of {A — XI) is dense for some A > 0. 

2. A is closable and its closure (also denoted by A) generates a contraction semigroup. 

We know that FI C Xi, and that FI is dense in Xi. By Proposition 12.81 we know that L 
generates a contraction semigroup on FI. As B is bounded on FI, we know that A (restricted to 
FI) generates a semigroup on FI. Thus it must be (see e.g. Theorem 1.5.3 in [T3]) that for A > 0 
large enough A — XI is invertible on H. Thus the range of A — AI contains H, and thus is dense 
in Xi. Because A is dissipative by ProDOsition l2.61 the Lumer-Phillips theorem then implies that 
A generates a contraction semigroup on Xi. □ 

By combining Proposition [2II] and Lemmas 12.41 and 12.91 along with Proposition 12.81 we imme¬ 
diately obtain the following: 

Theorem 2.10. The closure of L generates a semigroup on Xi uniformly bounded in time: 

We^^Waxi) < M forallt>0. 
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3 Polynomial Decay Estimates 


In this section we prove Theorem 11.11 namely that perturbations of equilibrium small in Xk will 
decay with a uniform, polynomial rate. We will first prove polynomial decay results for The 
following interpolation result is a modification of a theorem in [5], where it was originally used 
to study the convergence of travelling waves. 

Theorem 3.1. Let rj £ (0,1) and m, fc £ M with 0 < m < k. Let {S'(t)}t>o be a family of linear 
operators on Xi which for any t > 0 satisfies 

||S'(t)u-5'(t)u||xi <M||u-u||xi, ll'S'WwIlx, < 

where u, v are arbitrary elements of the appropriate spaces, M is a fixed positive constant and 
Xrj > 0. Then the operators S(t) necessarily are bounded from Xi^k to Xi+m and satisfy 

\\S{t)u\\xi+„, < C{1 + ||u||xi+fc for all u £ Xi+k and t > 0, 

where C depends on m, k, M and A,,. 

Proof. The proof is very similar to that found in [8], with modifications necessary, however, to 
handle the mass constraint and weighted norm on Xi. 

1. Consider : R x Xi —)• R defined by 


K{s,u)= inf (llw-t'IUi+e'*||u||xj- 

In interpolation theory [4] this is known as a modified K-functional. For fixed s, K(s,-) is a 
norm. Clearly K{s,u) is increasing in s and bounded above by ||M||jsfi- Furthermore, we claim 
that K is absolutely continuous in s. Indeed, if we define K(s,u) := iF(logs, u), then K[-,u) can 
be written an the infimum of affine functions, and thus must be concave. This readily implies 
that K{s,u) is absolutely continuous in s. 

We begin by proving upper and lower bounds on K. First, we get the lower bound 

OO OO 

K{s, u)>'^Qi inf (|ui - v\i + = 'Y] Qi\ui\{i A . (3.1) 

vGM. 

2=1 2=1 

Next, observe x A = x for all real x if and only if s > := —1 — log 77 . Thus for s > Sj^, 

OO 

K{s,u) < ||m||xi == X] ^ 6®+’'*) . 

i=l 

Suppose now that s < s^. Then l /?7 £ {cc : < a;} = [z_, z+] C (0, 00 ). Let j{s) be the least 

integer greater than or equal to z+, and define the sequence Vs{u) by 


i(u)i '.= < {Qil-) E Qkkuk 

k>j(s) 


for i < j{s), 
for i = j{s), 


0 


for i > j{s). 
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In particular note that QiiVs{u)i = 0, so Vs{u) € Xjj. Writing j = j{s), we then find 


K{s,u) < ||?/-?;s(?x)||^, + 

i>j 

/ pS+r}j \ ^ 3 — ^ 

< f 2 + j X] ■ 

' ^ ' i—j 2—1 

Now, , and i > j implies i = i A e®’*"’**. Furthermore, whenever 

1 < i < Z- we have e®’*'’*® < = Z- < llr] < ijr] = (i A e'*+''®)/? 7 . By these estimates we 

find that with C = max{2 + , 1/??} we have that for any s G R, 

00 

iF(s,M) < C^Q2k|(iAe*+®'®). (3.2) 

i^l 

2. In the next step, for r > 0 we set 


E 

*>i 




t-1 

E *3* 


i Qj 


,VJ 


Qjj 


oo 


QilUi 


hr{s) := 


{l-sY 


for s > 0 , 
for s < 0 , 


and define the norm 


111111*:= / K{s,u)hr{s) ds . 

Jr 


We claim this is equivalent to the norm in By (EH) and E3), it suffices to show there 

exist C -, C+ > 0 independent of i such that 


C'_(l + i)^+’'< {i A e^+’^YKis) ds < C+il + iy+^ fori>l. 


(3.3) 


To show this, we first bound the part of the integral over s G [0, oo), finding that 

pOO 

1< (iAe"+’’®)e-*ds<i< (l + i)^+^ (3.4) 

Jo 

For the part over s G (—00, 0], after changing variables twice via z = —s, a = z — rji, we have 

/ O poo 

{i A e*+’’®)(l - s)®-! ds<i (1 A e-^+®'®)(l + zf-^ dz 

-OO «/ 0 

pTji pOO 

= i / {l + zY~^dz + i / e“‘®(l + 771 + (t)®'“^ c?(T 

Jo Jo 


<Ci{l+r]iY <C(l + i) 


1+r 


This establishes the upper bound in (13.31) . 

To get the lower bound, choose so large that i > Irj implies rji — logi > ^rji. For i < /^ we 
have (1 + 1 )®'+^ < (1 + ^ 27 )’^''"^: hence we get the lower bound in (13.31) with C_ = (1 + 7,,)“^“®’ by 
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using (1X1) . For i > Irj, we find 


{i A - sy-Us = i (lA 


„ — 2 +i)i—log i 


)(1 + zY ^ dz 


>i z^-^dz>c{i+iy+r 

Jo 


Thus II • II* is equivalent to || • Hxi+r- 


3. Now, let Hr{t) := y°° hr(T)dT and k{t) := ^K{t,u). We claim that 

Hm{s+t)<cHkis)ii+tr-y 


for all s G K., and for t > 0. To prove the claim, we first note that 


Hm{s) 



for s > 0, 
for s < 0, 


and furthermore, for s < 0, we can find c, C > 0 so that 


c{l-sr <H^{s)<C{l-sr. (3.5) 

We then consider separate cases. First, if s > 0, 

Hmis + t)= < Ce-yi += CHk{s){l + 

Next suppose that s < 0 < s + A Then 

Hm{s + t) = < C(1 + s + !)-'= = < C(1 + ty'^Hkis), 

where we have used (13.51) . Finally, in the case that t < —s, we note that because m — k < 0, 

(1 - (s + t)r < (1 - sy < (1 - sYii + . 

In light of (13.51) this proves the claim. 

4. Next, we use the assumed bounds on our operators to estimate 

K{s,Sit)u)< inf (||5'(t)u-5'(t)u||xi+e®||5'(t)u||xj 

V£Xrj 

<M inf (||w-ii||xi+e'’"^''*||v||xj 
= MK{s — Xy, u). 

We remark that for u G we have that 0 < K{s,u) < HuUxi A e®||it||x„, and thus for u G 
we have that Hr(s)K{s, u) goes to zero as s ^ ±oo. Thus we may use integration by parts, and 
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our previous estimates, to obtain the following for any u £ 

\\S{t)u\\xi+^ <C [ K{s,S{t)u)hm{s)ds 

Jr 

< C / K{s — Xr)t,u)hm{s) ds 

Jr 

= C k{s^u)Hjn{s X-qt) ds 

Jr 

< C'(l+t)™"'=y' k{s,u)Hk{s) ds 

= C{l+t)^-^ [ K{s,u)hk{s)ds 
Jr 

= c(i+ tr->^\\u\u < c{i+ 

Because is dense in Xi+fc, we have the desired inequality. This completes the proof. □ 

We will apply the previous theorem to the semigroup generated by L. We first state a 
proposition, which will be an immediate consequence of Lemma 14.61 in the following section. 

Proposition 3.2. The operator L generates a semigroup on the space Xi+fc, for any k > 0. 

Proof. Applying Lemma 14.61 when g = 0, along with the fact that B is a bounded perturbation 
gives the desired result. □ 

With these tool in hand we can establish the following linear decay estimates. 

Corollary 3.3. Provided 0 < m < k, the semigroup e^* generated by the operator L satisfies 

||e^‘u||xi+,„ < C(1 + ||m||xi+, for all u G Xi+k, 

where C depends on m and k, but not on u or t. 

Proof. This follows directly from Proposition 12.81 Corollary 12.101 Theorem 13.11 and Proposition 

IQ □ 

Our goal is to use the detailed decay rates in Corollary 13.31 along with Duhamel’s formula, 
to prove Theorem 11.11 We first prove that Duhamel’s formula is justified in the appropriate 
spaces. We begin by recalling a fact from Ball, Carr and Penrose Proof of Theorem 2.2 and 
Proposition 2.4). 

Proposition 3.4. Let (ci) be a solution to the Becker-Doring equations, and let (hi) be defined 
by (11.61) . Suppose that h{0) G Xi_|_fc, with k >0. Then ||h(t)||xi+fc < ||h(0)||xi4.fcC'e^* for some 
C and K independent of h. 


We now have the tools to justify Duhamel’s formula. 
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Lemma 3.5. Assume that {ci{t)) is a solution of the Becker-Doring equations and {hi(t)) is 
defined by (11.61) . and let m > 0. If h{0) € Xa+m then the following is satisfied (strongly) in 

—h = Lh + hiTh. (3.6) 

In particular, if h{0) £ then we have that the following is satisfied in Xi^m- 

h{t) = e^*h{0) + f e^^*~^^hi{s)Th{s) ds, (3.7) 

Jo 

where is the semigroup generated by L on Xi+m (see Provosition VJ .2\] . 

Proof. Because h(Q) £ by Proposition 13.41 and Lemma [2.31 we have that Lh + hiTh is 

bounded in X 2 +m on any finite interval. Because each hi is continuous by definition (II.2|) . it 
must be that Lh + hiTh is measurable in X 2 +m- We claim that in X 2 +m we have that 

h{t) = h{Q) + f Lh{s) + hi{s)Th{s) ds. (3.8) 

Jo 

Indeed, the right hand side of the equation is well-defined, and must match the coordinate-wise 
integrals from definition 11.21 This implies that h{t) is locally Lipschitz in X 2 +m- As (13.81) also 
holds in we thus have that h(t) must be differentiable in This implies (13.61) . 

Again by Proposition [3^ we know that hiTh £ L^((0, T); Xi+m). Corollary 4.2.2 in [T3] then 
implies (13.71) . □ 

In deriving uniform bounds we will need a specialized version of Gronwall’s inequality. 

Lemma 3.6. Let u{t) be a positive, continuous function on [0,oo). Suppose that u satisfies 

u{t) < C2{1 + t)~’'+ f Ci{l+t — s)~^u{s)ds. (3.9) 

Jo 

Furthermore, suppose that r > 1 and that Ci is small enough that 

Cl [ {l + t-s)-^l + s)-^ds<0{l+t)-^ (3.10) 

Jo 

for some 6 < 1 and for all t > 0. Then we must have that 

Proof. Let v(t) = u{t)(l + ty. Then we have that 

v{t) < C 2 + {I + ty f Ci(l -I- f - s)“’'(l -I- s)“’'u(s) ds. 

Jo 

This then readily implies that for any T > 0, 


lkllc(o,T) <^2-1- 6 (||u||c(o,t)- 
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Thus for alH > 0 


which establishes the desired result. 


v{t) < 


C2 

1 - 6 »’ 


□ 


Remark 3.7. We note that for any r > 1 we can find a Ci > 0 such that (13.101) is satisfied. 
This is because 

ft ft/2 

/ (1 + s)-''(l + t - s)-^ ds = 2 (1 + s)-’’(l + t - s)-^ ds 

or+l 

<—-{i + t)-r 

r — 1 

Thus if Cl < {r — 1)2“^’’+^^ then we have that (|3.10l) is satisfied. 

Remark 3.8. The dependence on the constant Ci is critical in the previous proof. Indeed, if 
/q°° 1 ^ 1(1 + s)~^ds > 1 then it is possible to show that for some u{t) = c > 0 the inequality (13.91) 
is satisfied. Thus decay estimates can only be obtained if Ci is sufficiently small. 

The last tool that we need is a local stability estimate. The proof of this estimate is somewhat 
involved, and we postpone it until the next section. 

Theorem 3.9. Let (ci) be a solution to the Beeker-Doring equations (see Deftnition \1.2\} . and let 
(hi) be determined by (11.61) . Assume that the model coejficients in (11.21) satisfy (I1.11|) - (I1.14I) . Fix 
k > 2. Then given any e > 0 there exists S > 0 such that if ||h(0)||xi+fc < ^ then ||h(t)||xi+fc < e 
for all t > 0. 

With these tools in hand we now prove the main result of this paper. 


Proof of Theorem \1.1[ By Lemma 13.51 we know that the equation 


h(t) = e^*h(0) + f hi{s)Th(s) ds 

Jo 


is satisfied in Xi+m, where e^^ is the semigroup generated by L. By Corollary 13.31 we can thus 
estimate 


\\h{t)\\x,^^<C(l+t)-^’^-"^^\\h{0)\\x,^,+CUk,m) f\\hi(s)Th(s)\\xA^+t-s)-^’^-^-^'>ds. 

Jo 

By Lemma 12.31 we know that T is bounded from to and thus 


\\h{t)\\x,^^ < C{l + t)-^>^-^'>\\h{0)\\x,^,+CL(k,m)Cr f \\hi{s)h(s)\\x,^,{l+t-s)-^>^-^-^Us. 

JO 
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It is then immediate that 

\\h{t)\\x,^^ < C(1 + ||h(0)||x,+, + Csup \\h{T)\\x,^, [\l + t - ds. 

T Jo 

We then use a crude bound to obtain 

r Jo 

By Lemma [3.91 for any e > 0 we can choose Sk^m small enough to guarantee that 

\\hmx,^^ < C(1 + +£ f\i + t- ds. 

Jo 

As fc > m + 2, by applying Lemma 13.61 (whose conditions will be satisfied for e small due to 
Remark I3.7L we then find that 

which is the desired result. □ 


4 Local Stability Bounds in Xi+k 

In this section our goal is to prove Theorem l3.9l The general strategy is to derive bounds on the 
evolution family U{t, s) generated by F{hi{t)) when hi is small. We first establish bounds in H 
directly using dissipation estimates. We then establish stability bounds on U{t, s) in Xi+fc by 
using the extension principle from ProDOsition l2.ll This then immediately implies Theorem l3.9l 

4.1 Non-linear stability in H 

The following lemma gives a local, non-linear stability estimate in the space H. 

Lemma 4.1. Suppose that g[t) G for I = [0,r) with T possibly infinite. Suppose 

furthermore that the model coefficients in satisfy (jl.ll|) - (jl.l4|) . Then there exist Sh and 

A > 0 such that if \g{t)\ < Sh then {F{g{t))}t(^i generates an evolution family Uh in H on the 
interval I with bound 


\\UHit,s)\\c(H) < e for 0 < s <t <T. 

In order to prove this lemma, we will use the following proposition from Pazy (Corollary 5.4.7 
and Theorem 5.4.8 in [II]). 

Proposition 4.2. Let X be a Banach space and let I = [0,T), with T = oo permitted. Suppose 
that, for any fixed t € I, C(t) is the generator of a semigroup {Sc{t)is)}s>o which satisfies 


ll-S'c(t)(s)ll£(JC) < e foralls>0 
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where A is independent oft. Also suppose that dom(C'(t)) = D is independent oft and that for 
all X & D we have that C{t)x is in X. Then the family of operators {C{t)}t^i generates an 

evolution family U on X which satisfies 

\\U{t, forO<s<t<T. 

Furthermore for xq & D we have that x{t) := [/(t, 0)xo is the unique solution of the non- 
autonomous Cauchy problem 

^x(t) = C{t)x{t), a:(0) = a:o. 

dt 

Given fixed N, we define T to be a diagonal operator given by 

(Th)i = -aihi, := QiQi + bi, (4.1) 


define S to be the operator 

{Sh'^i .— ]1|2>A^+1} “b 11{2>A^}■ 

and K := L — T — S. In the proof of Lemma HIT] we will use the following facts (see Proposition 
2.10 and Corollary 2.11 in [^). 

Proposition 4.3. Assuming (Il.lll) - (ll.l4p . the operator L given by (II.9p satisfies the following 
properties: 

1. L is self-adjoint in H, with dom(L) = dom(T) = 

2. For some Ac > 0 we have that {h,Lh)^ < —Ac||h|||^ for all h € dom(L). 

3. L = T S K, K is compact on F[, and for N large enough, S is symmetric and satisfies 

< 6(||T/i||i/ for all h € dom(T), where 9 <1. 

We now prove Lemma l4.II 

Proof of Lemma \4.1\ We first claim that the following spectral gap estimate holds as long as g 
is sufficiently small: For some Ah > 0, 

{F{g)h,h)j^ < —AnWhW^ for all h G dom(L). 

To prove this inequality, we recall (ll.8|l and use Proposition 14.31 to estimate 

{F{g)h, h)^ = ((1 - e)Lh, h)^ + e {Kh, h)j^ + {{gV + e{T + S))h, h)^ 

< -(1 - e)Ac\\h\\l + e\\Kh^H,H)\\hrH + ((^r + ^(T + S))h, h)^ . 
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We select e small enough that > £\\K\\i^(^h,h)- As S is T-bounded with T-bound 9 < 1 

we have that S is relatively bounded (with relative bound smaller than one) by Because 

S is symmetric, this then implies (see [12], p.292. Theorem 4.12) that 




h, h 


< 0 . 
H 


Thus we can estimate 


_ (1 - £)Xc 

2 

< (l-g)Ac 

2 

< (l-g)Ac 

2 


WHh + — 2 —^ + h,h^ 

^ _ Q 

||/l||^ H" ^ ^ Qi ( ^ ^ hi ^l) 


ll^llff + ^ Qi 

i=l 



IQiffl 

2 





ll^llff + ^ Qi 

i=l 



1-6 

2 


tTi + aiC\Qig\ 


h 


2 
i 1 


where we have used the assumptions p.l2l) and (11.151) and the fact that Qi^i is finite. By 
(HU) there exists a. Sh > 0 so that if < 5h then {aiC\Qig\ — £^-^<7i) < 0. Thus if \g\ < 6 h 
we deduce that 

(F(g)h, h)^ < =: -XhMI, 

which proves the claim. 

We observe, from the previous estimates, that indeed ||r/i||ff < CUT/iH/i. This implies that 
S + gT is relatively bounded by T with relative bound strictly less than one, as long as |g| < 5h, 
where perhaps we have made 5h smaller. As T is self-adjoint, by Theorem 1.3.2 in [lOj we have 
that F{g) generates an analytic semigroup in H. Furthermore, by the relative bound it is clear 
that dom(F( 5 )) = dom(T) = P{Qi(Ji). 

Now, as g{t) is it is clear that for € H we have that F{g{t)))v is in C^{I]H). We then 
directly apply Proposition 14.21 to obtain the desired result. □ 


4.2 Non-linear stability in Xi+fc 

The main goal of this subsection is to prove the following lemma. 

Lemma 4.4. Suppose that g(t) G (^^(/jK), for I = [0,T) with T possibly infinite. Suppose 
furthermore that the model coefficients in satisfy ()1.11D - (I1.14I) and that k > 0. Then there 
exists a 5k such that if | 5 (t)| < Sk then {F{g{t))}tei generates an evolution family CNi+j, (t, s) in 
Xij^k on the interval I with bound 


\\Uxk(.'t,s)\\c{Xi+k) ^ A/fc, 


where Mk is independent of s,t and the particular choice of g. 
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To prove this lemma we will use Proposition 12.11 in conjunction with the stability in H 
established in the previous subsection. These techniques should also be applicable in the spaces 
but for the sake of clarity we do not pursue the analysis here. 

To begin, we define the operator A{g) in weak form by 

OO OO 

(4.2) 

— QN-lQiaN-lhN{,4>N — 4>N-1 — (kl) , 

where is a constant, greater than fV^ + l, to be determined. We then define B{g) := F{g)—A{g). 
The next proposition establishes the dissipativity of A{g). 

Proposition 4.5. Under the assumptions of Lemma \4-4\ if N in (14.21) is chosen large 
enough, then there exists a constant Sk so that if \g\ < Sk then 

{sgn{h),A{g)h)^.^^^^^_^^ < 0 for all h G X 2 +k- 
Proof. With Wi = and using (fi = WiSgii{hi) in (14.2p . we compute, as in the proof of 

Proposition mni 

(sgn(/i),^(<?)/i)x*^,.Xi+, 

OO 

= ^ Qihi(^QiaiWi^i{sgTi{hi^i) - sgn{hi)) hiWi-i{sgn{hi-i) - sgn{hi))^ 

i=N 

OO 

+ X! - U)i) + bi{wi-i - Wi)) 

i=N 

OO 

+ sgn(/ii) ^ Qihi{bi - Qiai) 

i^N 

OO 

+ ^ ^ QihiQiai{wi^iSgn{hi^i) - WiSgii{hi) - sgn(/ii)) 

i=N 

=: El + i?2 + E^ + i?4. 

First we estimate E 2 , written as 

OO 

i=N 

By choosing N sufficiently large we can make the ratio arbitrarily close to 1. Thus we 

apply (11.1611 to find that 

OO 

£'2 < -C ^ Qi\hi\ai{wi+i - Wi). 

i=N 


ciiQi - h, 


Wj - Wj-I 
‘ Wi+l - Wt 
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We next calculate 

OO 

E3 < Qi\hi\{bi + Qiai). 

i=N 

Recalling (I1.13p . and using that Wi+i — Wi ^ 00 since /c > 0, we thus have, for N sufficiently 
large, 

00 

E2+ E3 < -C ^ Qi\hi\ai{wi+i - Wi). 

i=N 

Because hi{sgn{hi±i) — sgn(/i,;)) < 0, we infer Ei < 0. Thus, in the case g > 0 we estimate 

OO 

El E 4 < E 4 < \g\ ^ QiQiai\hi\{wi 4 .i - Wi 1 ) 
i=N 

OO 

<q 5 lE Qi\hi\ai{wi4.i - Wi). 

i^N 

For g < 0 we find that 

E’l + E’4 < ^ Qi\hi\Qiai{-2wi^i - giwi^-i + Wi)) 

i>N 

sgn{hi)^sgn{hi+i) 

OO 

i=N 

When |(;| < 1 we have that the first term in (14.31) is negative. This then readily implies that for 
N sufficiently large and for \g\ sufficiently small we have that 

OO 

(sgn(/i), < -C'^Qi{wi+i - w,)a,\hi\ < 0, 

i=N 

which completes the proof. □ 

The next step is to prove that {^(^(t))} indeed generates an evolution family. 

Lemma 4.6. Suppose that the assumptions of Lemma EZ] are satisfied. Suppose furthermore 
that 

|g(t)| < min{4, Sk+iJn}, (4.4) 

where Sk is given in Proposition If. 51 and Sh in Lemma \4. 1\ Then for N chosen as in Proposition 
\4-5\ the family {A{g(t))}t^i generates an evolution family Rxi+fc on the interval I = [0,T) in the 
space Xi+k, which for 0 < s < t < T satisfies 


s)||£(js:i+fc) < 1- 

To prove this lemma, we use the following proposition, which is a direct application of Theorem 
5.3.1 in [H]. 
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Proposition 4.7. Let I = [0,T), with T = oo permitted, and suppose that a family of linear 
operators {C(t)}t^i satisfies the following for all t G I. 

1. C{f) generates a contraction semigroup on Xi^k- 

2. C{t) generates a contraction semigroup on X 2 +k- 

3. C(t) is a bounded operator from X 2 +k to Xi^k, cind the map 1 1 —> C(t) is continuous from 
I to C{X 2 +k, Xi^k)- 

Then {C{t)}tei generates an evolution family satisfying s)||£(Xi+fc) < 1- 

We now prove Lemma 14.61 

Proof of Lemma \4.6] We claim that {A{g{t))}t^i satisfies the assumptions of Proposition 14.71 
By (14.41) and Proposition 14.51 we have that A{g{t)) is dissipative on Xi^k and X 2 +fe. For fixed 
t S /, by (14.41) . F{g{t)) generates a semigroup on H (as established in the proof of Lemma [4. II) . 
As B{g{t)) is a bounded operator on H, it then must be that A{g{t)) generates a semigroup 
on H. This then implies that for some large, positive real A we must have that the range of 
^( 5 ( 0 ) ~ ^ contains H. Thus the range of A(g(t)) — A is dense in Xi+k and X 2 +k- As in the 
proof of Lemma 12.91 this implies that A[g{f)) generates a semigroup on Xi^k and X 2 +k-, and 
thus the first two assumptions are satisfied. 

Next, as g{f) is and by (jl.l4p . the third assumption is necessarily satisfied. Thus we may 
apply Proposition [321 which proves the lemma. □ 

We remark that the previous lemma is independent of the choice of g(t). The next result 
follows from a computation as in Lemma 12.41 and we omit the proof. 

Lemma 4.8. Under the assumptions of Lemma \4..6\ the operator B{g{t)) is uniformly bounded 
from Xi to H, with a bound that depends only on 5k, and not on g or t. 

We are now in a position to prove Lemma 14.41 

Proof of Lemma \4.4\ In light of Lemmas 14.61 and 14.81 this follows from Proposition l2.ll □ 

Remark 4.9. We note that the bound Mk is not dependent on the particular function g{t), and 
only on its bound 6k- This is because of the independence on g{t) in the bounds obtained in 
lemmas 14.61 and 14.81 

The following is a direct application of Theorem 5.4.2 in El- 

Proposition 4.10. Suppose that, for some k > 0, {C{t)}t^i generates an evolution family U in 
Ai_|_fc on the interval I = [0, T). Furthermore, suppose that for some h G C{F, X 2 +k) Xi+k) 

we have that 

jh = amt) 

is satisfied in Xi^k- Then it must be that U{t,0)h{0) = h(t). 
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Now we finally give the proof of Theorem 13.91 

Proof of Theorem I j’.ffl Let be the uniform bound in the space Xk given in Lemma 14.41 Set 

. _ Qi min{4_2,4-1,4,4+1, 

2Mfc 

Now, let {hi) correspond to a solution of the Becker-Doring equations, with ||/i(0)|lxi+fc < S. By 
Lemma 15751 and as fc > 2 we know that h^ is C^. By Lemma HiTl we thus know that {F{hi{t))}t^j 
generates an evolution family U on and on the (non-empty) interval / such that 

\hi{t)\ < min{(5fe_2, (5fc, (5fc+i, (5//}. As fc > 2, by Lemma 15751 we know that the conditions of 

Proposition 14.101 are satisfied in Ai_|_(j._ 2 ), and thus U{t, O)fc(O) = h{t) for all t € I. 

The uniform bounds from Lemma [4.41 then imply that ||fc(t)||xi+fc < -^*11^(0)11x1+*, the 
interval I. Our choice of S immediately implies that I = [0, oo) and that ||fc(t)||xfc < £/2, which 
completes the proof. □ 
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